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Abstract: The leading relativistic and QED corrections to the ground 
state energy of the three-body system e~e + e~ are calculated numerically 
using a Hylleraas correlated basis set. The accuracy of the nonrelativistic 
variational ground state wave function is discussed with respect to the con- 
vergence properties at the increase of the basis dimension and to the variance 
of the energy expectation value. Recent progress in the numerical procedure 
used to calculate expectation values for products of various physical opera- 
tors is presented. It is shown that the nonrelativistic ground state energy 
can be calculated with an accuracy below the level width. The corrections to 
this energy include the lowest order Breit interaction, the vacuum polariza- 
tion potential, one and two photon exchange contributions, the annihilation 
interaction, and spin-spin contact terms. The relativistic effects and the 
residual interactions considered here decrease the one electron binding en- 
ergy from the nonrelativistic value of 0.012 005 070 232 980 10(3) a.u. to 
0.011 981 051 246(2) a.u.. 
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I. Introduction 

The positronium negative ion (Ps~) is the simplest system composed of 
three equal mass fermions, e~e + e~ , bounded only by electromagnetic inter- 
actions. Similar examples of three-body systems, bounded by increasingly 
complex interactions, are provided by three-quark systems such as the proton 
and neutron, and three-nucleon systems, as the 3 H, 3 He nuclei. 

In a nonrelativistic approach, accurate numerical approximations to the 
bound eigenstates of three quantum particles interacting by Coulomb forces 
can be obtained by using the Rayleigh-Ritz variational method. A suitable 
set of coordinates and basis states for the three-body problem was proposed 
by Hylleraas [1] during the early days of quantum mechanics, and it was used 
to calculate the ground state energy of the helium atom. With respect to this 
set, the matrix elements of various two-body operators can be expressed in 
analytical form [2], and extensive high-precision calculations become feasible 

[3] [4]. ' ' ^ 

The relativistic quantm many-body problem can be approached either 
from the field theory, or by using a Schrodinger equation with an "action at 
a distance" type Hamiltonian, defined by quantizing the classical relativistic 
system [5]. Though, a puzzling result in classical mechanics is the no inter- 
action theorem [6], which apparently rules out any instantaneous action-at- 
a-distance Hamiltonian. This theorem states that in a classical many-body 
system, the relativistic invariance of the equations of motion (the physical 
laws) is compatible with the "manifest relativistic invariance", of the world 
lines determined by these equations, only if there is no interaction between 
the particles. However, this strong result can be avoided if the interacting 
particles have a structure, as the condition of manifest invariance becomes 
ambiguous [7] [8]. 

The approach to the bound state problem based on field theory leads to 
a relativistically invariant Bethe-Salpeter equation [9], p. 196. In the case of 
two relativistic electrons, approximate Lorentz invariance to the first order 
is introduced by the Breit interaction, which can be seen as the quantum 
correspondent to the Darwin term in classical electromagnetism [10]. 

In the helium atom, the two electrons move in the Coulomb field created 
by a composite, heavy nucleus, which to a first approximation can be consid- 
ered as center of mass (CM). The case of Ps~ is different, because all three 
particles are elementary, have the same mass, and move to the same degree 
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with respect to the CM. 

The existence of a bound ground state in the e~e + e~ system was pre- 
dicted by Wheeler [11] and was observed by Mills [12] passing a positron 
beam through a thin carbon film in vacuum. The measured Ps~ — > (27) e~ 
decay rate A = 2.09(9) nsec" 1 [12] corresponds to a Ps~ lifetime of 0.478 
nsec, intermediate between that of para (singlet) Ps (0.125 nsec) and ortho 
(triplet) Ps (140 nsec) [13]. 

Accurate nonrelativistic numerical calculations for the ground-state prop- 
erties of Ps~ are presented in refs. [14] to [17]. The autoionization states 
have been studied in [18], while several low-lying resonances have been pre- 
dicted recently [19], by using a combination between the stochastic varia- 
tional method (SVM) with correlated Gaussians and the complex scaling 
method. 

The accuracy of the Ps - groundstate wave functions given by SVM in 
a Gaussian basis, was studied by comparison to the direct solution of the 
Schrodinger equation in [16]. It was shown that despite the fact that in SVM 
the convergence properties of the expectation values for most operators are 
better, the wave function is less accurate. 

In this work, the accuracy of the Ps~ nonrelativistic variational ground 
state is studied by using beside the convergence properties of the energy with 
basis size, also the variance of the Hamiltonian. The numerical procedure 
used to calculate matrix elements is presented in Sect. II. It is shown that in 
agreement with [16], the variance is larger than the accuracy resulting from 
convergence. Estimates of the relativistic correction terms and the leading 
QED corrections are presented in Sect. III. Tables containing the expecta- 
tion values of singular operators appearing in the correction terms, such as 
p A and delta functions are given in the Appendix. The main results and the 
concluding remarks are summarized in Sect. IV. 

II. The Nonrelativistic Quantum Three-Body Problem 

The nonrelativistic Hamiltonian of the three-body system e~e + e" (or 
e + e~e + ) is 

Ho = (-£v? - ivl - /Vx • V 2 ---- + —) / a.u. (1) 

2 2 ri r 2 r 12 

where / = /i/m, \i = m/2 is the reduced mass, Vj = d? { , fi = Ri/a^, 
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i = 1,2,3 denote the position vectors in the CM frame of the two electrons 
(i = 1,2), and of the positron (i = 3) in units, while r = \r\ — f 2 |, 
T\ = \fi — r 3 \ and r 2 = |r 2 — r 3 | are the relative distances. The space 
coordinates unit is = a //, where a = % 2 /{me 2 ) = 0.529 177 249(24)A 
is the Bohr radius. By this choice, the Hamiltonian is naturally expressed 
in reduced atomic units of energy /a.u. (= 13.605 698 1(40) eV = lRy if 
/ = 0.5), where 1 a.u. = e 2 /a = a 2 mc 2 is the atomic unit of energy and 
a = e 2 /(hc) = 0.007 297 353 08(32) is the fine structure constant. 

Approximate eigenfunctions of this Hamiltonian are obtained by using 
the variational method. The trial function is a finite linear combination 

*(ri,si;f 2 ,s 2 )= ]T [q l ±(l,2)<S> abchhLM (ri,?2) (2) 

a,b,c,h,h 

+Qab! (2, l)$afe C i 1 J2LAf(r 2 ,fi)]^Sm.(si, S 2 ) 

of iV fe basis elements $ a6c h i 2 LM^Sm s - The orbital component $ a6c ; 1 ; 2 lm is 
represented by the Hylleraas correlated wave function [1] 

<Ka bc M 2 LM(n, f 2 ) = r^rJ 2 e-^-^^(l, 2) (3) 

which is a product between a polynomial in all relative radial variables and 
the orbital angular momentum eigenstates 

3^(1,2)= £ C l ^ M Y hmi (j- x )Y hm2 (r-2) , (4) 

m\+m2=M 

where f\ — fj/rj, i = 1, 2, r — r/r are unit vectors. 
The spin function 

lis 1 1 

^Sm s = Yl C ^i%2m s \-IIl)\-ll2) (5) 

corresponds to singlet (S* = 0) or triplet (S = 1) configurations, when the or- 
bital part is symmetric (^^(l, 2) = q 1 ^ (2, 1)) or antisymmetric 2) = 
— ?a b c (2, 1)), respectively. The expansion coefficients (1, 2), and the non- 
linear parameters a, /3, have been determined previously [17] by using the 
variational equations 

(gjgojjr) 
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The kinetic energy operator for the electron 1 is — Vf/2, and its effective 
action in the Hylleraas model space is given by 



V?$ = \-—r\ — + 1 _ 1 (7) 

r\ dr\ dr\ r 2 dr dr r\ 

+ 2(r,-r a f,-f.)_g_ _ Aft . vr^|* 

r dr\dr r\r dr 

where ii = fi x pi and W{ = —%r\ x / x . A similar expression, obtained by 
permuting the indices 1 and 2, yields V|$ . 

The action of such operators on the Hylleraas basis states is complicated, 
and in the case of their product becomes tedious, involving hundreds of 
polynomial terms rfr^r , spherical harmonics, and singular delta functions. 
Therefore, in this work the matrix elements have been calculated by using 
a new procedure, based on the representation of the physical operators as 
linear combinations within a set £ = {E k , k = 1, 2, n} of n elementary op- 
erators. Although S is not a Lie algebra, this approach makes the numerical 
calculation more flexible, because the components Ej. can be programmed in- 
dividually, and they can be assembled as needed to form various complicated 
operators. By a suitable choice, the set S may account for several physical 
operators of interest. Some of the elementary operators used in the present 
calculation are presented in Appendix, Table I. 

The accuracy of the wave function improves with the dimension N b of 
the basis set. When N b increases, the expectation value of the Hamiltonian 
(H )(N b ) = (^\H \^) decreases, and in principle, at the limit — > oo the 
series (iJo)(iv 6 ) approaches the exact ground state energy. Thus, a measure 
of the accuracy is provided by the deviation from 0, the known limit value, 
of the "convergence speed" [(iJ )(jv£') ~ (-^o)(Af^)]/(^6 / — where N' b and 
N{J are two consecutive values of N b . A different measure of the accuracy 
is the variance a = ^(Hq)^ — (H ) 2 Nb y The variational ground state en- 
ergy E g = (if )(;v 6 ) [17] and the present results obtained for a 2 are given 
in Appendix, Table II. The effective value defined in the Appendix- (b) is 
Egff = -0.524 010 140 465 960 215 38(56)/ a.u.. Previous estimates of E g 
in Ps~ by using the correlation-function hyperspherical-harmonic method 
and the stochastic variational method are —0.524 010 139 / a.u., respec- 
tively -0.524 010 140 452 / a.u. [16]. The effective value obtained in [15] by 
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using the extrapolation formula E g (N b ) = Ef f + A/N£ is -0.524 010 140 
465 956(8) / a.u.. 



III. Relativistic and QED Corrections 

The quantum description of a relativistic charged fermion based on Dirac 
equation requires two spin-1/2 wave functions, and corresponding 
to the retarded and advanced waves, respectively. For the particle eigenstates 
with energy E ~ mc 2 , ~ K\ip + ), K = a-p/(2mc), and the normalization 
condition + = 1, can be written in terms of the large compo- 

nents as = 1, where |0) = (1 + K 2 /2)\ip + ) . The antiparticle eigenstates 
are related in principle to the solutions with E ~ —mc 2 , when become 
the large components. Therefore, in general, the normalized eigenstate of a 
nonrelativistic Schrodinger equation for a particle or antiparticle should be 
seen as an approximation for |0) ps (1 + K 2 /2)\%pi) , where \ipL) denotes the 
corresponding large component. By definition, if Hi\ipi) = El\iPl), then 
Hl,k\4>) — El\4>), where H L K = H L + [K 2 ,H L ]/2 includes both the nonrel- 
ativistic Hamiltonian and the lowest order correction terms. 

In the interacting e~e + e - system Hl k will be restricted to H + Hi + H 2 , 
where Ho is given by Eq. (1). The term 

(.Pi — \Pi\) takes into account the relativistic variation of the mass with 
velocity, and 

H, = ~hprM,V]] (9) 

1=1 

derives from the sum between the term 2(pV) -p)/ (2mc) 2 of Hi and the com- 
mutator [K 2 ,V]/2 = 2[p 2 ,V}/(Amc) 2 . Here V(R, R u R 2 ) = e 2 (l/R-l/R 1 - 
1/R 2 ) = Vf a.u. and R = a^r, R 1>2 = a^r 1>2 . 

The magnetic current- current interaction plus the retardation correction 
corresponding to the lowest-order Breit interaction are described by the ad- 
ditional term 

M 2 = -^^{R^lPi ■ P2 + r ■ (f • p ± )p 2 ] (10) 
-R{ 1 [p 1 ■ p 3 + ri • (n • pi)pz\ - R 2 X [P2 -P3 + r 2 - {h ■ p 2 )p 3 }} 
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such that the total Hamiltonian for Ps that will be considered here is H — 
H + H 1 + H 2 + M 2 . 

In a classical relativistic many-body system, the dynamical CM defined 
by the condition J2iPi — 0, is not necessarily the same as the geometrical 
centre of mass, located at Rcm — J2i m iRi/ J2i m i- 111 the present case, 
if if = Hq + Hi + H 2 + M 2 , only the dynamical CM is inertial, because 
[pi + I>2 + p*3, H] = 0, while [fi + r 2 + r 3 , if] is not a constant. The use 
of the Hylleraas basis ensures that the present calculation takes place in 
the dynamical CM frame, because by the choice of the coordinates — iV» = 
Pi/ (fame) and V3$ = (— Vi — V 2 )§ for any variational wavefunction <3>. 

The expectation value (Hi) = — (1/64) (Vf + V 2 + Vf)a 2 f a.u. can be 
calculated either directly, or by assuming that in the ground state (H O p ) = 
(O p H ) = Eg(O p ) for any operator O p , and using the equalities Vf + V 2 = 
2(# - V + /Vi • V 2 ), with Hq = H /(f&.u.), and 

Vt + V 4 2 = A(H -V + fVi- V 2 ) 2 -2VM , (11) 

= (-Vi - V 2 ) 4 = + VI + 4(Vi • V 2 ) 2 (12) 

+2V?V^ + 4(V? + V^)Vi • V 2 . 

Although formally the same, within a finite basis the two expressions give 
slightly different results ((Vf), (Vf)s), presented in the Appendix, Table III. 
In the numerical estimates we have used only (Vf)#, because of its rapid 
convergence and higher acuracy in the effective value. 

The term H 2 contains the singular operators AiV = — 47r[5(r) — <5(ri)], 
A 2 V = -4n[6(r) - 6(f 2 )}, and A 3 V" = 47r[5(ri) + S(f 2 )}, which yield 

(H 2 ) = a 2 it(5(ri) + S(f 2 ) - 5(r))f a.u. (13) 

Previous estimates of (S(Ri)) (= a M 3 (5(r 1 ))) in Ps - by using the correlation 
- function hyperspherical - harmonic method and the stochastic variational 
method are 0.020 733 14(6)a 3 , respectively 0.020 731 048 976 a 3 [16]. The 
same methods give for (5(R)) the values 0.000 170 997(2)a 3 and 0.000 171 
112 600 741 a 3 , respectively [16]. The results of the present calculation, in 
the same units (a^ 3 ), are listed in Appendix, Table IV as a function of the 
dimension iVj, of the basis set. 

The expectation values which appear in the calculus of (M 2 ), obtained 
when N b = 324 are 

Uee = (r _1 Vi • V 2 ) = -0.008 267 646 67 , 
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Vee = ( r - l r ■ (r • Vi)V 2 ) = 0.019 610 925 35 

and for i = 1, 2 

Uep = (r-^Vi • V 3 ) = 1.535 434 049 31 

Vep = (r^n ■ (n ■ V 3 )Vi) = -0.555 009 821 912 . 

In terms of these variables, (M 2 ) = 0.5a 2 wf 3 a.u. with w = u ee + v ee — 
2u ep — 2v ep = —1.949 505 176 125. For the 2528-dimensional basis set 
w = —1.949 505 250 368. The average of the last three consecutive values, 
obtained for N b = 1990, 2276 and 2528 gives the effective matrix element 
w eff = -1.949 505 250 368(1). 

The effective sum of the spin independent relativistic correction terms 
(H^, (H 2 ) and (M 2 ) is -0.145 476 184 397(8)a 2 / a.u. which decrease the 
Ps~ ground state energy to 

E e /J = (H} eff = Ef f - 0.145 476 184 397(8)a 2 /a.u. . (14) 

The same calculations yield for the corrected ground state energy E®» in 
neutral positronium = — (0.5 + 5a 2 /32)/ a.u.. However, for this relativis- 
tic two-body system the finite mass corrections to the energy provided by 
the one-body Dirac equation can be obtained exactly up to the order a 2 Ry 
by using the formula [20] 

iW) = 7^-1 -{(??- I) 2 ] /a.u. (15) 

where r] = l/yjl + (Za) 2 /(n - u) 2 , v = j + 1/2 - yj(j + 1/2) 2 - (Za) 2 . The 
expansion £(1,1/2,1) ~ —(0.5 + 5a 2 /32)/ a.u. reproduces E®*, showing that 
the relativistic corrections given by the expectation value of H are reliable. 

Within QED the constituents of the three body-system e~e + e - cease to 
be "elementary", because they are subject not only to the mutual two-body 
Coulomb-Breit interaction, but are also coupled to the vacuum fluctuations 

— * 

of the electromagnetic field A [21]. The interaction terms accounting for 
this coupling are represented by an infinite series of increasingly complicated 
Feynman diagrams with closed photon lines. However, the complexity is 
increased recursively, by taking into account at each order three basic pro- 
cesses, represented by the anomalous magnetic moment (vertex) corrections, 
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electron self-mass and vacuum polarization diagrams. 

Although formally complicated, the main effect of the coupling to the 
field degrees of freedom is simply a change in the charge and mass param- 
eters e and m of the theory. This contribution has already been taken into 
account, because it is included in the measured values of e and m used to 
define the atomic unit of energy. Though, the QED corrections in the inter- 
acting three body-system e~e + e~ are not the same as for the free particles, 
and the differences still need to be considered. 

The vacuum polarization properties have been studied first by Heisen- 
berg [22] and Uehling [23], showing that a given charge density p(R) induces 
a polarization charge Sp(R) = — (a/157r) Aq V 2 p(i?), where Ao = h/mc is the 
Compton wavelength of the electron. The induced charge leads to deviations 
from the standard Coulomb interaction. Thus, the vacuum behaves as an 
inhomogeneous dielectric, in which the mutual potential energy between two 
point-like charges Z\ and Z 2 is [23] 

V{R) = ^-[l--RU{R)} , (16) 

by U(R) denoting the Uehling potential. This potential is singular at it! = 
0, falls of exponentially for R > 0, and satisfies the integral condition 
/ d 3 RU(R) = — 47tAq/15. Therefore, it can be well approximated by a delta 
function, U(R) = — Air (X^/ 15)5 (R). In the case of Ps~, the correction term 
introduced by this potential is 

(H vp ) = ^a 3 (5(f) - Sin) - 5(f 2 ))f a.u. (17) 

By using the effective values given in Appendix, Table IV, the contribution 
of the vacuum polarization to the Ps~ ground state energy is [H vp ) e ** = 
-0.022 024 212 934 6{7)a 3 / a.u.. It is important to remark that this value 
takes into account the positron recoil (the "mass polarization" term) because 
the wave functions are obtained by minimizing the full nonrelativistic Hamil- 
tonian. In neutral positronium, (5(Ri))p s = l/ijia^) = l/(87rao), and the 
vacuum polarization correction is — l/(157r) a 3 f a.u.. 

As it was shown early by the Lamb shift measurements [24], the main 
QED correction appears however from the coupling to the vacuum fluctua- 
tions of the field rather than from the vacuum polarization ( [21] p. 59). For 
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a free electron the ground state energy is given by its rest mass m = uib + Sm, 
consisting of the uncoupled value mj, and the positive renormalization term 
5m = (3amb/2n)\n(A/mb) due to the electromagnetic self-energy, where A 
is a large (formally infinite) cutoff mass. 

Similarly, the coupling to the field modes also affects the intrinsic excita- 
tions of a many-body system. In a bounded iV-particle system, the shift AE n 
in the energy E n = (n\H \n) of the level \n) due to the exchange of a trans- 
verse photon can be obtained by using the time-independent second-order 
perturbation expression 

AE n = -(n,O f \H c (A)——^ H c (A)\n,O f ) . (18) 

Here H C (A) = X)£Li hi(A) is the sum over all particles of the one-body cou- 
pling terms hi(A) = -e^ • A^, di « Pi/(m,ic), 

^ = / 5 E Wl*-** (19) 

» A — 1)2 

is the quantized transverse vector potential of the photon (e^ • k — 0, = 1), 
-"'a = / d 3 k J2x=i,2 ^ lc ^ a k\ a k\ is f ne f ree fi e ld Hamiltonian, and |0j) denotes 
the photon vacuum. This shift has the form AE n = Z)£Li X£ + Si<j Y^, 
where 

Xf = -<n, OfMA)— hl (A)\n,0 f ) (20) 

and 

^ = -2^,0/1^——^ -h^faOf)] . (21) 

-no + tlA — J^n 

It is important to remark that the interaction with the vacuum field 
fluctuations may affect not only the intrinsic dynamics, but also the center 
of mass. In a classical two-body system coupled to the field, H c can be 
written in terms of the canonical pairs (r, p^) = (Ri — R 2 , \ip\jm\ — npijm'i) 
and (Rcm, Pcm) = {[xR\ / m 2 — nR 2 1 Pi +P2) of intrinsic and respectively, 
center of mass variables, as 

Pfi ■ [— A(r 2 ) ~ — A( Rl )\ —Pcm • [eiA (Rl ) + e 2 A [R2) } . 

Ul 2 ' TOi 171x1712 
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This expression shows that in a neutral two-body system (such as Ps) the 
center of mass energy is not affected by the field only if A( Rl ^ = A( R2 ), or 
when the size of the system is negligible compared to the photon wavelength 
(dipole approximation). 

In a quantum iV-body system it is convenient to take advantage of the 
finite size effects by writing A as the incoherent sum of long and short 
wavelength components, Ai and Ag, obtained by decomposing / d 3 k as 
f\k\<k L d 3 k + J'k L< \k\ < k M y w h ere an d ku are cutoff parameters. Each domain 
brings its own contribution to the matrix elements, which can be similarly 
decomposed as 

Xf = Xf" + Xf" , Y^ = Y^ + Y^ n . (22) 

At the end of the calculation k L should disappear, while k M — > oo. 

If H consists of the kinetic energy term plus a local potential V, then a 
non-relativistic calculation within the dipole approximation yields 

X i n = -ir^Y-^hckMfM) + (n\[p t ; [p t , V]]\n) ln^ - 253] , (23) 
6nm z c z k R 

where k R = Rm/^c, Rm is a dimensional constant with units of energy, and 
Bfi are the diagonal elements of the matrix [B™] defined by 

B ij = E(^n ~ £ m )Re((ra|fllm) • (mjp») In |K ~ . (24) 

The first term depends only on the kinetic energy, and it can be written 
as — 8rriL(n\pf |n)/(2m 2 ), 5m L = 4mr e A;i/(37r), where r e = ah/(mc) denotes 
the classical radius of the electron. It contributes also to the energy of a 
free particle (V = 0) and has the structure of a first-order perturbation shift 
induced by a variation bm^ of the nonrelativistic mass. Thus, such terms 
can be taken into account simply by a redefinition of the cutoff mass A. 

A relativistic calculation of the one-body QED correction arising from the 
exchange of a transverse hard photon at a Coulomb vertex [25], [21] p. 177, 
yields 

x -^s?^y>^ • < 25 > 

(5/6 = 11/24 + 3/8) such that 

A? - -|>Wl»> + + f)«M0. + m . (26) 
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The quantity Y^ n can be expressed as 
^ = -^<fi-S>.- S ^ ? ^ i «K-,B,V]]>.ln^-2 fl5 }.(27) 

In the case of Ps~ there are three terms Xf, one for each electron (i = 
1,2) and one for the positron (i = 3), and three terms Y£ n , i < j. The 
contribution to AE n arising from the terms linear in SrriL of X and Y L is 
5oE n = —5mL(n\(pi + p 2 — p3) 2 |n)/(2m 2 ). In the dynamical CM frame, this 
energy shift can be accounted, for example, by an effective variation 3SmL 
in the total mass of the electron-electron pair and SrriL in the mass of the 
positron, or by a variation of ASmi, in only one of them. 

The definition of the Bethe logarithm f3 n = 2B% 3 / (\p 3 , -[p 3 , H ]]) n , and 
the identity m 2 J2 i , j e i e j B n (i,j)/(m i m j ) = e 2 J B n (3,3)(l + m/m 3 ) 2 (valid if 
mi = rri2 = m and (n| Y^iPi\n') = for any n,n f ), show that Si^f + 
E l<i y;f* = 5 £„ + b x E n + <5 2L £ n , where 

r\h mc Fi _ 

4 A = W?^^" + (ln 2ft^ + 6> S^"" (28) 

= ^-[-4,a„{i 3 (r 1 )+i 3 (r 2 )>„+2(hi-^ + |){i 3 (r 1 )+* 3 (r 2 )-i 3 (r)>„]/ 3 aaL 
3 z/zkr o 

and 5 2 LE n = Z)j<i(^ / ij')n is given by the expectation value of the potential 

WS(fcL) = ^ln^5 3 (ry)/ s a.u.. (29) 

The term Y$ n due to the exchange of a short wavelength (hard) transverse 
photon between different particles will be decomposed as Y^ n = Y^ 2n + Y^ 3n , 
according to the expansion 1/(H + H A — E n ) pa l/H A — (H — E n )/(H A ) 2 . 
The contribution from 1/ H A is 

Y** = -2( n MH&)Hl%{As)\"^f) . (30) 

In the limit ki — > 0, ku — > oo, the integral over k in this matrix element can 
be evaluated by using the identity 




■ B — k ■ Ak ■ B) = —(A ■ B + r ■ Ar ■ B) , 



12 



showing that the sum J2i<j YE 2n becomes the two-body term (M2) n of order 
a 2 Ry, already taken into account. Thus, the only new contribution is the 
next-order term 



Y** = 2Re[(n,0 f \h i (A s ) I ^-^h j (A s )\n,0 f )] (31) 
which is the expectation value of the two-body potential 

Wg(k L ,k M ) = ^f(k L: k M:rij ) +AnS 3 (r ij )\n^]f a.u. (32) 



Here f(k L ,k M ,r) = 2[j (k M r) + j 2 (k M r) - jo(k L r) - j 2 (k L r)]/(3r 3 ) is the 
function introduced by Araki [26], written in terms of the spherical Bessel 
functions jo 5 j2- When fc L — > 0, k M — > 00, /(0, 00, r) = — 2/(3r 3 ), but the 
logarithmic factor in the second term of Wf, is divergent at both limits. 
However, the divergence in k^ is cancelled by the low-energy term, and the 
sum W i3 \k M ) = W£(k L ) + Wg(k L: k M ), 



W l3 {k M ) = ^mk M ,r l3 ) + %lln^]/ 3 a.u. , (33) 
37r I r Xi k R 

is independent of ki- The divergent factor containing kM contributes only 
when \n) is an S state, but in this case the expectation value (l/r 3 ) n is also 
logarithmically divergent. It is however possible to define a limit for the sum 
of these infinite terms in the sense of the principal value. Let 

^/ x 9(r — a) 5(r) , a u 

D(a, r) — — — - — In — 34 

r 6 r l a 

be a distribution depending on the positive radius parameter a = r]/kM, 
where 77 is a positive scale factor. Because r 2 d a D(a,r) = [5(r) — 5(r — a)]/a, 
when kM — > 00 the expectation value (D(a,r)) n is finite. In terms of this 
distribution we can define the principal value 

V[h(0, k M , r) + ^ In ^] | fcM _oo = 4vr5 3 (r) In - lun a ^ Q D(a, r) . 

(35) 

4 

The choice of a scale factor r/ = es -7 , where 7 is the Euler's constant, yields 
the formula used by Araki [26] 



2a 



3 



(Wij) n = - — {Qij + An(5\r l3 )) n [\na»k R - -]}/ 3 a.u. , (36) 



3tt 
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where Qfj = lim a ^o(D(a,r) + 4ix r )5 3 (rij)) n . In Ps this yields for the effec- 
tive two-body contribution 8 2 E n = S 2L E n + Ekj^ln = Ei<j(Wy)n the 
expression 

2(T 3 

S 2 E n = - — [Q n 12 + Q n 13 + Q n 23 (37) 

+4n(\na,k R - ^)(n\5 3 (n) + 5 3 (r 2 ) + 5 3 (r)\n)]f 3 a.u. . 

Summarizing these results, the effective QED contribution of order a 3 to 
the energy level E n of Ps~ due to the exchange of a transverse photon is 
Si p E n = 5±E n + 5 2 E n . This sum is independent of the arbitrary energy unit 
Rjh, as it should, but for the purpose of numerical calculations we choose 
R-M — /Ry- With this choice, a^k^ = a/2, and mc/(2hk R ) = l/(fa 2 ). 

The corresponding terms for positronium can be obtained from the ex- 
pressions given above simply by neglecting all the expectation values con- 
taining the variables r 2 and r, involving the second electron. For the Ps 
ground state Q^ 3 = —4 In 2, while f3g S is the same as the Bethe logarithm for 
hydrogen (3% = 2.984 128 555 765 497 611(4), each Bethe logarithm being 
calculated using the corresponding reduced Rydberg constant [27]. For the 
Ps~ ground state the numerical values of Q{ 2 and Q{ 3 used in the present 
estimates are listed in Table V, while g = 3.005 030(2) [27] (including the 
finite mass correction). 

To the same order we should consider also the double photon exchange 
term (including the Coulomb part) S 2p E n [26] [28] 

a 3 

5 2p E n = ~^[Qi2 + Q13 + Q23 ( 38 ) 
4 1 3 

-47r(ln/a- - In 2 + — )(n\5 3 (n) + 5 3 (r 2 ) + 5 3 (r)\n)]f 3 a.u. , 
3 6 

and the energyy shift associated with the two-photon decay. 

In general, any coupling which makes the levels unstable produces a com- 
plex energy shift A c E n = S c E n — where 5 c E n is a correction to the 
level centroid, A" = T™/% is the decay rate, and c denotes the decay channel. 
Neutral positronium normally decays by spontaneous e + e~ annihilation into 
two photons if the total spin S ep = 0, (S ep = s e + s p ), and in three photons 
if S ep = 1 [29]. The corresponding decay rates are such that T 37 ~ aT 27 , 



14 



and the first correction arises from two photon annihilation. In this channel 
5 2y E n /T^ = -(1 - hi2)/7r [30], where 

T n 2l = 2na 3 ((2 - S 2 ep )5 3 ( ri )) n f a.u. . (39) 

For the Ps ground state (S ep = 0) this yields a decay rate \(p s ,2~f) = a 3 Ry/h = 
8.04nsec _1 , in good agreement with the experimental result 7.99(ll)nsec _1 
[31]- 

In the Ps ground state the electron spins are coupled to 0, and the two- 
photon annihilation can take place between the positron and any of the two 
electrons. The total rate depends on (Sf 3 + S 2 3 ) g = 3, and can be expressed 
in the form 

r (Ps - >27) = 27ra 3 (5 3 (r 1 )>/ 3 a.u. . (40) 

The effective ground state expectation value (S(Ri/a )) given in the Ap- 
pendix, Table IV, yields A(p s - >27 ) = 2.092 797(1) nsec" 1 , in good agreement 
with the previous estimates [32] and the experimental result 2.09(9) nsec -1 
[12]. The corresponding level shift is S 2l E g = —(1 — ln2)r(p s - 27 )/vr. 

Summarizing the results of these calculations, the effective ground state 
expectation values of the nonrelativistic Hamiltonian and the first relativistic 
and QED corrections for Ps and Ps~ are collected in the following table 



E[Rj] Ps Ps" E Ps -E Ps - 

E g -1/2 -0.524 010 140 465 960 215 38 0.024 010 140 465 960 215 38 

±0.56 x 10~ 18 ±0.56 x 10~ 18 

(#i)/a 2 -5/32 -0.161 254 673 938 50(6) 0.005 004 673 938 50(6) 

(H 2 )/a 2 1/4 0.259 466 645 837(8) -0.009 466 645 837(8) 

{M 2 )/a 2 -1/4 -0.243 688 156 296 0(1) -0.006 311 843 704 0(1) 

(H vp )/a 3 — 1/(15tt) -0.022 024 212 934 6(7) 0.000 803 553 855 7(7) 

S lp E g /a 3 2.766 873 00(3) 3.006 491 9(9) -0.239 618 9(9) 

S 2p E g /a 3 -0.585 335 778(7) -0.510 831 605(7) -0.074 504 17(1) 

8 2l E g /a 3 -(1 - hi2)/7r -0.025 448 161 055(1) -0.072 226 124 976(1) 



In Ps the nonrelativistic one electron binding energy 0.024 010 140 465 
960 215 38(56) Ry is practically the same as the one determined in [15], 
and close to the older estimate of 0.024 010 113 Ry [32]. The effect of the 
correction terms is to decrease slightly this energy to 

B' = [0.024 010 140 465 960 215 38(56) - 0.010 773 815 602(8)« 2 (41) 
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-0.385 545 7(9)« 3 ] Ry = 0.024 009 416 924 85(6) Ry . 

In positronium the observed hyperfine splitting of 1.160 963(9)a 2 Ry be- 
tween the otherwise degenerate ground state components corresponds to an 
additional spin-spin contact term [30] 

4 rv 1 -, 9fi n 

8?K = 2™ 2 (5 3 (r 1 ){- S - e - S - p (l--) + -5 e 2 p [l-( y +ln4)-]» n f a.n. , (42) 

shifting the energy of the singlet by SfE® = — 2a 2 (l — «/27r)/ 3 a.u. = 
—0.265 947 576(23) x 10~ 4 Ry. In Ps~ the electron-electron spin-spin de- 
pendent energy shift is [26] 

5fE g = -ya 2 (l + A a )<Sk • s 2 5 3 (r))f a.u. , (43) 

while the two electron-positron spin-spin contact terms contribute by 

5^E g = ^(5^)01, + Sf 3 )[l - (f + In 4)^]) f a.u. . (44) 

These formulas yield an additional shift of the Ps~ ground state energy 
5 s E g = 8fE g + 5f p E g = 0.207 196 744(18) x 10~ 4 Ry, and a contribution to 
its binding energy of 8f E° g - 5 s E g = -0.473 144 32(3) x 10^ 4 Ry. Including 
the spin-spin contact terms, the binding energy becomes 

B" = B' + 5 e s p E° g - 5 s E g = 0.023 962 102 492(3)Ry . (45) 



IV. Summary and Conclusions 

The calculation of the relativistic and QED corrections to the energy lev- 
els of a quantum three-body system represents a challenging problem the 
modern theory. Difficulties appear both at conceptual and computational 
levels, as there is no satisfactory relativistic many-body quantum theory, 
and the nonrelativistic problem is not integrable. 

A quantum three-body system thoroughly investigated since the early 
days of quantum mechanics is the helium atom. In this system a major sim- 
plification occurrs, because the reduced electron mass /x is smaller than the 
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mass of the positive charge by a factor 1.3707 x 10~ 4 , and to a first approx- 
imation the relative motion of the nucleus in the center of mass frame can 
be neglected. The relativistic invariance is partly restored by the Breit inter- 
action, and highly accurate nonrelativistic wave functions can be obtained 
numerically, from variational calculations. Within this framework, a pertur- 
bative treatment of the relativistic and QED corrections gives energy levels 
in remarkable agreement with experiment [3] [4]. 

The same procedure was applied in this work to the negative positronium 
ion. However, by contrast to helium, all three particles have equall mass, and 
a perturbative treatment of the positron motion becomes inappropriate. 

The accuracy of the nonrelativistic energy and variational ground state 
wave function was discussed in Sect. II. The effective value of E g estimated 
here is -0.262 005 070 232 980 107 69(28) a.u., the same as in [15], [17] and 
close within 10~ 8 to estimates obtained by other methods [16]. The variance 
of the Hamiltonian for the largest (2528-dimensional) basis set is 2.78 x 10 -8 
a.u., smaller than the level width T = Ti\ = 5.06 x 10 -8 a.u. due to the 
(2^)e~ decay. 

The calculation of the first relativistic and QED corrections has been 
presented in Sect. III. Formal expressions of these correction terms have 
been known for a long time, but the mass polarization term, which cannot 
be neglected for Ps~, increases dramatically the complexity of the numerical 
calculations. In this work the calculations have been performed by repre- 
senting the action of the physical operators on the Hylleraas basis states in 
terms of a suitable set of elementary operators (Appendix, Table I). The 
numerical values obtained for some of the most important matrix elements 
are summarized in the Tables II-V of the Appendix. It was found that 
the spin independent relativistic terms contribute to the Ps~ ground state 
energy by -0.072 738 092 198(4)a 2 a.u. and the lowest order QED correc- 
tions by 1.224 094 00(44)o; 3 a.u.. These terms decrease the ground state 
energy to E e JJ = -0.262 008 467 959 9(4) a.u.. Both contributions decrease 
slightly also the one electron binding energy, from the nonrelativistic value 
0.012 005 070 232 980 10(3) a.u. to 0.012 004 708 462 43(3) a.u.. A much 
larger contribution appears however from the spin-dependent contact terms, 
which raise the ground state energy to E e /J = -0.261 998 108 122(1) a.u., 
and further decrease the binding energy to 0.011 981 051 246(2) a.u.. The 
calculated decay rate by two photon emission is 2.092 797(1) nsec -1 , close to 
the previous theoretical results and to the measured value [12]. 
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The binding and ground state energies are shifted also by relativistic cor- 
rections of order a 3 Ry containing spin-orbit magnetic interactions. These 
contributions are also important, and should be included in the further at- 
tempts to improve the accuracy of the present estimates. 

Appendix 

(a) Elementary operators 

In terms of the elementary operators defined in Table I, the nonrelativistic 
Hamiltonian contains linear combinations such as 

~ \^l)<& = (E 2 + E 3 -E 4 + E 5 + E 6 - 0.5E 9 

-0.5-Eio - 0.5-En - 0.5E 12 )<5> 
for the kinetic energy part, = — £13$ for the potential energy and 

Vi • V 2 $ = (E l + E 2 + E 3 -E A + E 5 + E 6 + E 7 + E 8 )§ 

for the recoil term. The action on the basis states of a linear combination 
Op = J2k fpkEk, has the general form 

t=i 

where x denotes the whole set of indices of $, q px is a factor determined by x 
and the action of Ej, on the radial functions, while A t denotes the remaining 
angular operator. In this sum the same state A t Q Xt (or A t A t r$ x , in the 
case of a product) may appear several times with different scalar factors q. 
Therefore, in general it is possible to reduce the number of terms from N p 
to N r < N p by partial summations, before the effective calculation of the 
matrix elements. The reduction increases the speed of the numerical calcula- 
tion, because iV r can be significantly smaller than iVp. For example, N p : N r 
is 20 : 18 in the case of Vi • V 2 , 400 : 235 for (Vi • V 2 ) 2 , while for H and 
Hi is 31 : 24 and 961 : 427, respectively. 
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In the computer program the action of consists of both arithmetic and 
symbolic operations. The arithmetic operations are determined by its action 
on the radial function r^r^r^e'^ 1 " 131 " 2 - The number n k of distinct polyno- 
mial terms generated by this action is given in the second column of Table 
I. The symbolic operation corresponds to the action of on the spherical 
harmonics in 3^lm, and is coded by a character denoting angular operators 
such as f\ • f 2 or Vf • V^f. 



(b) Error estimates 



The series of numerical values presented in Tables II- V appear to be con- 
vergent, but for comparison with experiment, it is useful to provide also a 
single effective value, representing the expected result of the present calcu- 
lation when N b — > oo. The procedure adopted here to define this value 
depends on the manner of convergence. In the case of a sequence {f n } 
convergent as an alternating series, the effective value / e // ± a/, given in 
the last row, was defined as the arithmetic average of its last three con- 
secutive terms, by f eff = (f nx + f ny + /„J/3, n x < n y < n z , and a) = 
[(fn x ~ feff) 2 + (fn y - feff) 2 + (/„, - /e//) 2 ]/3. If {fn} approaches the limit 
by monotonous increase or decrease, then we can assume that the series can 
be extended to infinity by the function F(n) = f e ff + Ae~ in . The matching 
equations F(n x ) = f Ux , F(n y ) = f Hy , F(n z ) = f Uz between F(n) and the last 
three calculated numerical values yield the parameter f e ff in the form [17] 

f - f 

/_ p | J J Tlx 

eff = Jriy + 



R-l 



Here R = e^ n y~ nx) is the solution of the equation R-l = /yi-^K-^VK-™-)], 
where 

fn v fn 



j ) ■ J n x 

tu, — 



fn z fn y 

The error is assumed to be 

I fn v fn 



°7 = \friy ~ feff\ = 



R-l 



If n y — n x = n z — n y , then R = R y . When n is simply Nb, then n y — n x = 286 
is larger, but close to n z — n y = 252, and R = R y still provides a reasonable 
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estimate. 



(c) Expectation values of p A 

The expectation values (Vf) = (V 2 ) for electrons in the Ps~ ground state 
can be calculated numerically either directly, as 

(V?) = A((E 2 - 0.5E 4 + E % - 0.5£ 9 - 0.5£ n ) 2 ) 

or as (Vj} E = 0.5(Vf + Vfj E = 0.5((V? + V 2 2 ) 2 ) E - (V 2 V|), where the first 
term is expressed in the form 

0.5<(V 2 + V 2 2 ) 2 )e = 2((E g - V + /Vi • V 2 ) 2 ) 

= 2[E 2 g - 2E g (V - fV 1 • V 2 > + ((V - /Vx • V 2 ) 2 )] 

(Eg = Eg/ fa.u.) by assuming that the variational ground state is practically 
eigenstate of Hq. Although formally the same at the limit JV& — > oo, the nu- 
merical values obtained for (Vf) and (Vf) E at finite AT 6 are slightly different. 
These estimates are given as a function of the basis size N b in the first two 
columns of Table III. The third column contains the relativistic correction 
term for the positron, given by 

(V|) B = ((Vi + V 2 ) 4 > £ = 2(Vi) E + 2(V 2 V 2 > - 8((E g - V)Vi • V 2 ) . 

(d) QED corrections of order a 3 in the limit m 3 — > oo 

When m 3 — > oo the vacuum polarization term and the contribution of the 
electron-electron spin dependent contact interaction 8 e s e E n \ a z = — (20/3)(s*i • 
s 2 8 3 (r)) n a 3 f 3 a.u. remain the same, but = 0, Y£ = 0, and Si p E n becomes 

KK = « 3 \\(^ ^ + I - M{i 3 (n) + <s 3 (r 2 )}„ 

-I< ln m~ R -\ + '"« 3 (''»» - 1^ 

The two-photon contribution reduces to 

5 2p E™ = a 3 [-^ + 2(\n fa - - In 2 + -)(5 3 (r)) n ]/ 3 a.u. , 
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and the total correction 5E™ = 5 lp E™ + 5 2p E™ + 5fE n \ a z + (H vp ) n is 

4 1 Q 1 4 1 64 7 

5iC = a 3 [-(--ln/a 2 -^)(5 3 (r 1 )+5 3 (r 2 )) n +( y ^fa+— )(5 3 (r)) n --Q" 2 ]/ 3 a.u. . 



Table captions 

Table I. Elementary operators {Ek, k = 1, 15} used in the calculation of the 
expectation values, denotes the number of distinct polynomial terms gen- 
erated by the action of E^ on a basis state. 

Table II. The ground state expectation values E g = (H ), (Hq) and a 2 = 
(Hq) — E 2 as a function of the basis dimension iVj,. 

Table III. Ground state expectation values of the singular differential opera- 
tors (Vf)s, (Vf) and (V|)_e as a function of the basis dimension N b . 
Table IV. Ground state expectation values of the singular Dirac distributions 
(5(Ri)) and (5(H)) as a function of the basis dimension N b . 
Table V. The ground state expectation values Qf 2 and Q{ 3 defined by Qfj = 
lim a ^o(^(^y — a )/ r fj + 47r[7 + ln(a/a^)]<5 3 (rjj)), as a function of the basis 
dimension N b . 
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Table I 



k 




E h 


1 


4 


fi • ToB 2 
i 1 ' 2 L/ rir2 


2 


4 


r^frofi ■ fo — ri)<9 2 

' V z 1 ' Z ' 1/ y\T 


3 


4 


4 


1 


d 2 r + 2r- 1 <9 r 


5 


1 


ri(r2r)~ 1 fi • V^f<9 r 


6 


1 


r 2 (rir) _1 f2 • Vf <9 r 


7 


4 


r 2 X ri ■ V\d ri + rf l r 2 


8 


1 


(nr^^Vf-Vr 
^ + 2rr 1 9 n 


9 


3 


10 


3 


d 2 r2 +2r^d r2 


11 


1 


r7 2 Vf • Vf 


12 


1 


r 2 - 2 vr • vr 


13 


3 


r 1 1 + r 2 1 — r 1 


14 


1 


-/i(/i + 1K 2 


15 


1 


-/ 2 (/ 2 + l)r 2 - 2 
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Table II 



N b 








E g [Rj] 












(#o>[V] 




a 2 x 10 14 [Ry 2 ] 


324 


-0 


.524 


010 


140 413 


399 


000 


28 





.274 


586 632 449 


596 


519 352.588 780 183 9 


411 


-0 


.524 


010 


140 455 


551 


566 


88 





,274 


586 628 565 


868 


126 562.113 711 645 6 


512 


-0 


.524 


010 


140 464 


139 


040 


54 





,274 


586 627 626 


769 


31 752.227 055 472 07 


630 


-0 


.524 


010 


140 465 


665 


621 


87 





,274 


586 627 375 


932 


6 508.536 436 039 192 


764 


-0 


.524 


010 


140 465 


918 


375 


12 


0. 


.274 


586 627 323 


102 


1 199.051 083 460 867 


918 


-0 


.524 


010 


140 465 


954 


391 


13 


0. 


.274 


586 627 313 


704 


255.472 732 778 647 


1089 


-0 


.524 


010 


140 465 


959 


038 


66 


0. 


.274 


586 627 311 


885 


73.086 762 209 023 06 


1283 


-0 


.524 


010 


140 465 


960 


002 


45 


0. 


.274 


586 627 311 


421 


26.589 554 852 763 07 


1495 


-0 


.524 


010 


140 465 


960 


160 


85 


0. 


.274 


586 627 311 


222 


6.672 153 792 304 994 


1733 


-0. 


.524 


010 


140 465 


960 


203 


19 


0. 


.274 


586 627 311 


175 


1.965 916 264 831 472 


1990 


-0 


.524 


010 


140 465 


960 


212 


96 


0. 


.274 


586 627 311 


165 


0.965 692 558 621 058 


2276 


-0 


.524 


010 


140 465 


960 


214 


82 


0. 


,274 


586 627 311 


160 


0.465 897 417 244 748 


2528 


-0 


.524 


010 


140 465 


960 


215 


25 





,274 


586 627 311 


158 


0.308 152 352 372 668 


eff 


-0.524 010 140 465 960 215 4(6) 


0.274 586 627 311 156(4) 


0.23(23) 













Table 


III 






N b 










(Vf) 


<Vf>£ 




324 


2. 


532 


451 


004 442 6 


2.532 


445 719 29 


5.255 396 862 


891 


411 


2. 


532 


451 


050 420 6 


2.532 


451 697 56 


5.255 397 122 


254 


512 


2, 


532 


451 


056 877 


2.532 


450 741 84 


5.255 397 117 


353 


630 


2. 


532 


451 


009 132 


2.532 


449 964 61 


5.255 397 051 


034 


764 


2, 


532 


451 


018 719 1 


2.532 


450 992 21 


5.255 397 086 


467 


918 


2. 


532 


451 


022 453 6 


2.532 


451 056 52 


5.255 397 094 


127 


1089 


2. 


532 


451 


021 529 7 


2.532 


451 023 43 


5.255 397 091 


672 


1283 


2. 


532 


451 


020 589 3 


2.532 


451 019 49 


5.255 397 091 


024 


1495 


2. 


532 


451 


020 595 


2.532 


451 022 17 


5.255 397 090 


993 


1733 


2. 


532 


451 


020 587 2 


2.532 


451 020 24 


5.255 397 090 


958 


1990 


2. 


532 


451 


020 559 2 


2.532 


451 019 92 


5.255 397 090 


940 


2276 


2. 


532 


451 


020 559 6 


2.532 


451 020 43 


5.255 397 090 


949 


2528 


2. 


532 


451 


020 560 


2.532 


451 020 42 


5.255 397 090 


945 


eff 


2.532 451 020 559 6(3) 


2.532 451 020 2(2) 


5.255 397 090 945(4) 
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Table IV 



N b 




(S(Ri))[ao 


3 ] 




(5(R))[a 


3 ] 


324 





,020 733 174 


230 


2 


0.000 171 000 


000 8 


411 





,020 733 203 


838 


1 


0.000 170 999 


383 2 


512 





,020 733 199 


804 


5 


0.000 170 999 


967 2 


630 





,020 733 193 


292 


2 


0.000 170 997 


306 7 


764 





,020 733 197 


986 


7 


0.000 170 996 


885 4 


918 


0. 


.020 733 198 


238 


9 


0.000 170 996 


811 


1089 





.020 733 198 


094 


3 


0.000 170 996 


832 4 


1283 


0. 


.020 733 197 


999 


5 


0.000 170 996 


756 


1495 


0. 


.020 733 198 


024 


3 


0.000 170 996 


767 3 


1733 





.020 733 198 


007 


4 


0.000 170 996 


760 1 


1990 





.020 733 198 


003 


4 


0.000 170 996 


757 7 


2276 


0. 


.020 733 198 


005 


3 


0.000 170 996 


757 1 


2528 


0. 


.020 733 198 


005 





0.000 170 996 


756 8 


eff 


0.020 733 198 004 6 


(8) 


0.000 170 996 756 7(4) 



Table V 



N b 


QU% 3 ] 






Q% 


K 3 } 




324 


0.095 757 780 


75 


-2. 


,776 


563 


295 


411 


0.095 757 975 


79 


-2. 


,776 


588 


343 


512 


0.095 757 804 


27 


-2. 


,776 


583 


829 


630 


0.095 758 749 


78 


-2. 


,776 


578 


687 


764 


0.095 758 904 


03 


-2. 


,776 


582 


810 


918 


0.095 758 930 


40 


-2. 


,776 


582 


894 


1089 


0.095 758 918 


40 


-2. 


,776 


582 


776 


1283 


0.095 758 949 


78 


-2, 


,776 


582 


694 


1495 


0.095 758 944 


76 


-2, 


,776 


582 


722 


1733 


0.095 758 947 


86 


-2, 


,776 


582 


703 


1990 


0.095 758 949 


04 


-2, 


,776 


582 


700 


2276 


0.095 758 949 


31 


-2, 


776 


582 


702 



eff 0.095 758 949 4(3) -2.776 582 702(1) 
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